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All-loop order critical exponents for massless scalar field theory 
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Abstract 

We compute analytically the all-loop level critical exponents for a massless thermal Lorentz- 
violating 0{N) self-interacting scalar field theory. For that, we evaluate, firstly explicitly up 
to next-to-leading order and later in a proof by induction up to any loop level, the respective 
/3-function and anomalous dimensions in a theory renormalized in the massless BPHZ method, 
where a reduced set of Feynman diagrams to be calculated is needed. We investigate the effect 
of the Lorentz violation in the outcome for the critical exponents and present the corresponding 
mathematical explanation and physical interpretation. 
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I. INTRODUCTION 


It has been shown that symmetry has been a very efficient and elegant gnide in the con- 
strnction process of physical theories. The Standard Model (SM) of high energy elementary 
particles is based on the U{1) 0 SU{2) 0 SU{3) internal symmetries. Another symmetry, 
a space-time one, to regard in the making of the SM is the Lorentz symmetry. The former 
are associated to the internal symmetry space of the helds and the latter to the space where 
the helds are embedded. In recent years there were many proposals to extend the SM to 
another realms. One of these attempts cnlminated in the so called Standard Model Exten¬ 
sion (SME), where one of the symmetries, the Lorentz symmetry, is broken [l-4|. Symmetry 
properties are of vital importance also in low energy physics, specihcally in phase transitions 
and critical phenomena. In his pioneering work, Wilson evaluated the loop quantum correc¬ 
tions to the classical Landau critical ex^nents in a massless Lorentz Invariant (LI) O(A^) 
self-interacting scalar held theory |^, by applying a theory invariant under Lorentz 

symmetry and transformations of the Euclidean Wdimensional orthogonal group. The crit¬ 
ical exponents can be computed in a held theory approach for this statistical mechanic’s 
problem. This approach is named thermal held theory jb]. In this approach the Lagrangian 
density plays the role of the energy of the system in the canonical ensemble and the par¬ 
tition function is similar to the generating functional of n-point correlation functions. For 
magnetic systems, the magnetization is proportional to the mean value of a huctuating held. 
As it is known, for rather distinct physical systems as a huid and a ferromagnet, the critical 
exponents are the same, showing their universal character. They can be written as functions 
of their dimension d, N and symmetry of some A^-component order parameter, and if the 
interactions of their constituents are of short- or long-range type. A few works exploring the 
ence of the critical exponents on the parameters with an easier experimental access 


depenc 

n 

as d [7|, l8[ anc 
parameter 


N jg 


11| as well as on a less accessible one, i.e. the symmetry of the order 


13| were investigated. We say that two or more systems belong to the same 
universality class when their critical behaviors are described by an identical set of critical 
exponents. In empirical terms, the O(A^) universality class is a generalization of the real 
models with short-range interactions: Ising {N = 1), XY [N = 2), Heisenberg [N = 3), 
self-avoiding random walk {N = 0) and spherical {N —>■ oo) [l^. As the mass of the held 


corresponds to the difference between an arbitrary temperature T and the critical one T^, 
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massless and massive theories describe critical and noncritical theories, respectively. The 
thermal fluctuations responsible by the corrections to the mean held values for the critical 
exponents are now related to the loop quantum contributions to the anomalous dimensions 
of the theory. The exponents are then extracted from the scaling properties of the renormal¬ 
ized n-point correlation functions near (massive theory) or at (massless theory) the critical 
point. So, equal results for them must be obtained if they are computed at massless or 
massive theories renormalized at the same or at different renormalization schemes. In fact, 
the critical exponents were computed in a few distinct theories renormalized at different 


renormalization schemes 


15| . Two of them are associated to the held and composite held 


anomalous dimensions, namely rj and respectively. The remaining critical exponents 7 , a, 
(3 and 5 can be evaluated by using four independent scaling relations among them. A critical 
theory is described by a LI Lagrangian density composed of a LI kinetic operator 
and a self-interacting term. The composite held correlation functions are generated by 
a term like where 0 ^ = (p‘^{y) = 0 (?/) 0 (|/) is evaluated at the same point of space-time. 
Terms with higher powers like (9^09^0)^, 0®, 0®... are not included because 


it is known that they give negligi 


are called irrelevant operators 


Die corrections to the critical behavior of the system and 
16l-ll8|. Operators with odd powers as 0^, 0® are ruled out 
by internal symmetry principles. As the magnetization of the system is invariant under a 
simultaneous change of all spin orientations, the Lagrangian density must be invariant un¬ 
der 0 —>■ —(f). An obvious extension to this theory it to consider its Lorentz-violating (LV) 
version and to investigate the consequences of this symmetry breaking mechanism on the 
outcome for the critical exponents. Such a theory is attained if we add to the LI Lagrangian 
density, a LV kinetic term of the form K^^d^cf/d^cf) 19|, where the dimensionless constant 
coefficients are not invariant under Lorentz transformations, i.e. 7 ^ Af^A'fKpcr. If 
\Kpiy\ 1, we have a slight violation of the referred symmetry. These coefficients are sym¬ 
metric {Kpy = K^p) and equal for all N components of the field and preserve the 0{N) 
symmetry of the internal symmetry space of field. Thus the introduction of the LV kinetic 
operator makes a connection between the breaking of an internal symmetry of the order 
parameter and a possible modification of the respective universality class, the 0{N) studied 
here. The achievement of this work is investigating the consequences of Lorentz symmetry 
breaking on the evaluation of the critical exponents. 

This work is organized as follows: In section [TTl we compute analytically, by explicit 
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calculations up to next-to-leading order, the critical exponents for critical thermal Lorentz- 
violating O(A^) self-interacting scalar held theory. In section HTTI the results of the section 
mi are used as inspiration to prove by induction the corresponding critical exponents up to 
any loop level. We hnalize this work in section HVl with our conclusions. 

II. NEXT-TO-LEADING LEVEL CRITICAL EXPONENTS 

The critical behavior of the system is described by the bare d-dimensional Lagrangian 
density 




+ + 2^b4> 


B- 


( 1 ) 


The quantities, (j)B, ^b and ts are the unrenormalized held, coupling constant and composite 
held source, respectively. In a thermal held theory approach, for getting a renormalized the¬ 
ory it is sufficient to renormalize the primitively bare T^^ T^^ and T^’^^ n-point correlation 
functions or IPI vertex parts 2nl-l22l|. When written perturbatively up to next-to-leading 
order (as the tadpole in a massless theory is null, the two-point function up to next-to- 
leading order must be expanded up to three-loop level), the unrenormalized T^^ and 
T^’^^ vertex parts have the expressions 


p(2) _ 

^ B — 


-1 




( 2 ) 


- X - iO -|- 2 perm. — - XXX + 2 perm. 

- + 5 perm., 



The theory is renormalized multiplicatively 






ryTlj^ ryl 


Qji 


( 3 ) 

( 4 ) 

( 5 ) 


where A = and g are the renormalized dimensional and dimensionless coupling constant, 
respectively and p is an arbitrary momentum scale parameter. In the BPHZ method, the 
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divergences are subtracted by counterterms order by order in perturbation theory. Thus we 
obtain the renormarlized theory 

L = ( 6 ) 

through the eq. [T]with 


^=zy''yB, 


and 





t = 



The renormalization constants in eq. [6] are given by 


i=l 


m 

( 8 ) 

(9) 


( 10 ) 


= i + ( 11 ) 

i=l 


Z^2 - l + 'Y^ 


2 = 1 


( 12 ) 


The c^, Cg and c ^2 coefficients are the i-th loop order renormalization constants for the 
held, renormalized coupling constant and composite held, respectively. The renormalization 
constants up to next-to-leading order are given by [23 1 


Zct>{g, 


.-- 1 ) 






-z he- 


(13) 


(XIX + 2 perm.) Sq + tc (XIX3<+ 2 perm.) S^xx 


Z,(<,,£-‘)-l+—^ 

+ 5 perm^ 


(14) 
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(15) 

where S q is the symmetry factor for the referred diagram and so on when some N- 
component held is considered. Now we proceed to show how to attain a full solution of 
the problem by evaluating not all diagrams, but just a few of them. The minimal set of 
diagrams for that task are 


€^--/ 


d^qi 1 1 

{27iY {27iY ql + K^^q^q'( + K^^q^q)^ 


1 

^ (gi + g2 + py + K^yqi + q2 + PY{qi + q2 + PY 


(16) 


= _\3 [ d‘^q2 d^qs 1 _1_ 1 

"O' J {27rY {27rY {‘^ttY qI + K^^q'^q^ ql + K^^q^ql ql + K^^q^q^ 

1 
X 


(gi + g 2 + Py + K^Yqi + q2 + PYYi + Q 2 + Py 

1 

(di + ds + py + K^Yqi + gs + PYYi + Qs + PY ’ 


O =>~Y 

where Pi + P 2 = P, 


d^q 


X5 


= -A 


{2nY g2 + Kf.^.qf^q'' (g + Py + K^Yq + PYY + PY ’ 

3 f d^qi d^q2 1 1 

{2nY {2nY ql + K^^qU'i {P - qiY + K^^P - qiY{P - qiY 

1 


X- 


(g^ + Kf^tyqf^q’^ + m?Y (Y + 


(17) 


(18) 


ql + K^^ql^ql (gi - q2 + PsY + K^Yqi - q2 + PsYYi “ d2 + P3Y' 

Writing the free propagator as an expansion in the small coefficients 

1 If K^^q^^q^ n{n + 1) K^^Kp^qf^q’^qPq^ 

1 — ^ H-;- , ^ - — -h ... 


(19) 


g2 _|_ ^2 


2 ! 


(g^ + m 


2\2 


( 20 ) 
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and by applying the formulas in the Appendix and dimensional regularization and e- 


expansion tools in e = 4 — d 


24|-|26| we obtain 




2d 2 r 


8e 


1 + -e - 2eLs{P) 


„2p2 

n" + ^AVir(P). 


( 21 ) 




6 e^ 


1 + -e — ZeL^i^P) 


^2 

2e 




-K,uL^/{P), 


( 22 ) 


O 

1 




i-y-i.L(P) 


n - -eK^,L"''(P) 




(23) 


where 


X5 


■ 2e2 


1- -e-eL(P) 


- eK^,L^'^{P) 


+-eK,,Kp„[L^^-{P)6P^ + L^'^P^{P)] \ , 


(24) 


n = 1 - + ..., 

Z o 

^ + 6PPS‘''^ + and 6^^ is the Kronecker delta symbol, 


L{P)= / dxln 
Jo 


x{l — x)P‘^ 

p 2 


(25) 


(26) 




Iq X {1 — x)P2 


(27) 


L‘i{P) = / da;(l — x) In 


a;(l — a;)p 2 

p 2 


(28) 


ir(P)= 


/o a;(l — a;)P 2 


(29) 


L-'<»(P)= f\jHi-=^rp“P‘'p^p<’ 
Jo 


(30) 


[a;(l — a;)P 2]2 

The diagrams O and XII evaluated perturbatively in the small parameters up 

to 0{K^) and the and ones up to 0{K). The computation of the later couple of 
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diagrams up to O^K"^) would be very tedious 


27 


. The others diagrams can be written 


in terms of the above as 




3^3p2 

16e^ 


1 + je - 2eL,{P) 


^4- 


8 e 


K,.L‘t:(P), 


(31) 


XXX = -A=/ 

f d‘^q2 


d'^qi 


{2tiY qf + Kfi^qUi {qi + Pf + K^Yqi + PY{qi + pY 

^ ^ 4 (or 


(2vr)‘^ ql + K^^q^q!^ (ga + PY + K^Yq^ + P)>^{q2 + PY " 

'‘‘® 111 - e - eUP)] n" - iK^L‘“'(P) + hA;„A-,,l2L'“'(P)i'’" + L'“"’"(P)| \, 


(32) 


0-0 

3/i"g3 




2 e 2 


1 - -e--eL{P) 


n2 - -eK^^L^^YP) + -eK,,KpY^L^YPY^^ + L^^'^^YP)] , 


(33) 


= A 


d'^q 


{2'kY g2 + Kf.^qf^q’^ {q + QY + K^Yq + QYY + QY ^ 


O 


P^Q 


1 - if - \eUQ) 


n - i€AVL'‘''(Q) + Y<i‘P<m\L‘“'(Q)S'° + r'“''”(Q)] I". ( 34 ) 


where Q = Qi + Q 2 , 

-S' = 

f d'^q2 


r d'^qi 1 1 

(27r)'^ qf + K^^q^q’( (gi + QY + K^Yqi + Q)^{qi + 
1 1 


(2vr )‘^ql + Kn^q2q2 {q2 + QY + Kfj.p{q2 + Q)^{q2 + Q)^ 


tXXX 


P^Q 


|1 - e - eL(Q)] n" - eK^,Li“'(Q) + -iK^[2Li“'(Q)6f^ + L‘“'>’-’(Q)] , 


( 35 ) 




























_ 2 /" 1 1 

j {2nY {2nY qj + (Q - q^y + K^^Q - qiY{Q - qiY 

1 1 


Q 2 + ^)^iy(l2l2 ~ (I 2 + QsY + K^^{qi — q2 + QY}^{,qi — q2 + Qs^ 

■-ML- 

~i{ n^-eK,.L‘“'{Q) + YK,.K^l2L'“'{Q)S>“’ + L'^^‘"{P)]\, 






1 - - ^e-t(O) n" - ieAV„L'“'(< 3 ) + YKp^KM‘“'{Q)6 '" + i'“''"( 0 )l} , 



I -g^-Mq 


1^1 l-i^-i^LW) tf-i£AVL""(Q) +j£/f,./fp,l2L"''(QM'" + L'‘“'"(Q)]}, 


The coefficients and c ^2 used above are the one-loop corrections to their respective renor¬ 
malization constants. They arise from the condition that the theory be hnite at one-loop 


* =-a'94 = --k ) 0 < , 


As the two-point function is hnite at one-loop order, = 0. Thus, we obtain for the 0{N) 
theory the renormalization constants 


Zd, — 1 


(A + 2)(7V + 8) A-ldny. 

1446 12966^ \ 4 / 


iV + 8„ , r(Af + 8)2 5/V + 221 . j 


Z,-l + -^Ilg + 


nv, 


9 




(43) 


— 1H—-—n^f + 
6e 


{N + 2){N + 5) N + 2 


36e2 


24e 


2^2 


^^9 


Observe that the momentum-dependent integrals in eqs. 12611301 disappeared in the hnal 
expressions for the renormalization constants. These constants depend on only through 
the LV factor If as a power-law. This fact will be essential for the developments of the 
next section. The renormalization constants were obtained for arbitrary external momenta, 
showing the generality of the method. The scaling properties of the IPI vertex parts are 
determined by the renormalization group equations for these correlation functions 

where 


Pig) = h 


9/i’ 


(45) 


7<a(^) = h 


dlnZ^ 

dfi 


(46) 


dlnZfi 

Mu) = 


(47) 


Now we evaluate the /3-function and held and composite held anomalous dimensions. They 
have the following expressions 


Pig) = 9 



iV + 8 
6 


lig- 


3N + 14 
12 



(48) 


i^ig) 


= 2 _ (iV + 2)(iV + 8) 

72 ^ 1728 

, , N + 2^ r 5 ^ \ 
'Ic^^ig) = —^^9 (^1 - Y2^9J ■ 


3^3 




(49) 

(50) 


The eqs. 14811501 are hnite functions of e, as promised by the renormalization program. The 
hxed points of the theory are computed from the eigenvalue condition Pig*) = 0. The trivial 
one or the one known as Gaussian is g* = 0 and is responsible for the classical Landau values 
for the critical exponents. For these exponents, the thermal huctuations at all length scales 
are not taken into account. These huctuations are considered when g hows to the nontrivial 
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fixed point. Its value computed from the parentheses in eq. 
given by 


up to the order of interest is 


9 = 


6 e 


(iv + 8)n 


1 + e 


3(3A^+ 14) 


(A^ + 8)2 


9 


*( 0 ) 


n 


(51) 


where is the LI nontrivial hxed point |20|. Thus by applying the relations between the 
critical exponents t], v and the held and composite held anomalous dimensions rj = 

= 2 — 'jrf>^{g*), respectively, the LV 11 factor cancels out and we obtain 

{N + 2)e‘^ r r6(3iV + 14) T 


9 = 


2{N + Sy 


1 + e 


(A^ + 8)2 




u = 


(N + 2)e (N + 2){N^ + 23N + 60)e2 _ 

4(Ar + 8)^ 8(lV + 8)=> “ 


(52) 


(53) 


where the exponents and are the corresponding LI ones 291]. In next section we 
evaluate the same critical exponents for any loop level. 


III. ALL-LOOP ORDER CRITICAL EXPONENTS 

The results of the earlier section show that after the cancelling of the coefficients 
through the momentum dependent integrals, the renormalization constants, /9-function, held 
and composite held anomalous dimensions and the hxed point depend on only due to 11 
in a power-law form. We can see from the eqs. ITTMdl 148115(1] and EH that the LV theory can 
be attained from the LI one by the substitution g^^^ —)■ Il^f in their loop quantum corrections, 
one power of 11 for each loop order, where g^^'' is the LI dimensionless renormalized coupling 
constant. In fact, it is known that the massless LV theory can be written in terms of 
the LI one with the effective scaled coupling constant = Il^f by applying coordinates 

19|. Thus, for any loop order, we 


transformations ideas [15[ inspired in the massive case 
have 


f^{9) = 9 




n=2 


7(9) = ET’nV. 


n=2 


lA'j) = 


"9” 


(54) 


(55) 


(56) 


n=l 
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where /3n\ and 7</,2 are the LI nth-loop corrections to the referred fnnctions. The Gans- 
sian fixed point is nnll as earlier in the next-to-leading level approximation. The nontrivial 
hxed point, evalnated from parentheses in eq. [511 is then = Il^f*, where g*^^^ is the LI 
all-loop order nontrivial hxed point. The LV 11 factor cancels out in the computation of the 
LV critical exponents at any loop level and asserts that the LV critical exponents are equal 
to their LI counterparts. This hnal result, that the critical exponents are the same both 
the LV and LI theories can be completely understood only if we evoke some ideas from the 
research branch of phase transitions and critical phenomena. This research branch is a very 
traditional and rich one in providing a large amount of experimental and phenomenological 
evidence that completely distinct system as a huid and ferromagnet present an identical crit¬ 
ical behavior. In other words, they share the same set of critical exponents. This happens 
when they have in common their dimension d, N and symmetry of some iV-component order 
parameter if the interactions are of short-range type. As the symmetry breaking mechanism 
does not occur in the internal symmetry space of the field and consequently of the order 
parameter but in the space-time where the field is embedded, the critical exponents must 
not change. In terms of the renormalization group language, for asserting the influence of 
an operator on the critical behavior of a system, it is not suffice to including it just based on 
naive dimensional analysis. Its final effect on the critical exponents must be checked after 
the renormalization process. 


IV. CONCLUSIONS 

In this work, we have analytically evaluated the all-loop order critical exponents for a 
critical thermal Lorentz-violating 0{N) self-interacting scalar field theory in the mass¬ 
less BPHZ method. For that, they have been computed in a proof by induction up to any 
loop level. The induction procedure was guided by an explicit calculation up to next-to- 
leading order. It was shown that the massles BPHZ method is simpler that its massive 
version in attaining the same goal, where a minimal set of diagrams was needed. We in¬ 
vestigated the effect of the Lorentz violation in the outcome for the critical exponents and 
showed that the critical exponents are identical to their LI counterparts. We furnished both 
mathematical explanation and physical interpretation of the results as a requirement for a 
full understanding of them. The mathematical explanation if that the LV coefficients 


12 


can be absorbed in an effective LI theory with a scaled conpling constant by applying coordi¬ 
nates transformations techniqnes. The physical interpretation is: as the Lorentz symmetry 
braking mechanism does not occnr in the internal symmetry space of the held bnt in the 


in which the held is embedded 


30[ |. the O(A^) symmetry is preserved. Thns the critical 


exponents must not change. This work sheds light, through a formal treatment, on the 
study of a possible change of the critical exponents under a symmetry breaking of the order 
parameter. 


Appendix A: Integral formnlas in d-dimensional Euclidean momentum space 

Considering Sd = Sd/i^irY = [2‘^“^7r‘^/^r((i/2)]“^ where Sd = lV{dl2) is the unit 
d-dimensional sphere area, we have 


r d^q _1__ ^ i r(d/2) T{a-d/2) 

J {27r)^ {q'^ + 2pq + '^2 r(a) (M^ — ’ 


r ir(d/ 2 ) p^r(a —d/ 2 ) 

J {27i)^ [q'^ + 2pq + M'^Y ^^2 r(a) ’ 


/ 


ddq qiiqu ir(d/2) 

(27r)'^ -|-2pg-|-*^2 r(a) 

'1^ r(a-l-d/2) ^ T{a-d/2) 

2 (M2 - p2)a-l-d/2 ^ ^ (jVfa _ p2)a-d/2 


1 


(A2) 


(A3) 


/ 


(Yq q^q^qP 

(27r)'^ (g2 -)- 2pq + M2)" 


o i r(d/2) 

'^2 r(a) 

r(a - 1 - d/2) 

(M2 - p2)a-l-d/2 


- 1 - p^^p’^pP 


r(Q; — d/2) 

(M2 - p‘^Y-dl2 ’ 


(A4) 


d'^g 


qpq’^qpq^ 


{2tiY {(f + 2.pq -I- M‘^Y 


= Sd.r: 


ir(d/2) 


X 


2 r(a) 

Y +0 0 J('j\g-2 _p2)a-2-d/2 




r(a-d/ 2 ) 


+pPp'^pPp" 


(M 2 - p^y-dji 


(AS) 
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